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Aims. The gas-phase abundance of methanol in dark quiescent cores in the interstellar medium cannot be explained by gas-phase chemistry. In 
fact, the only possible synthesis of this species appears to be production on the surfaces of dust grains followed by desorption into the gas. Yet, 
evaporation is inefficient for heavy molecules such as methanol at the typical temperature of 10 K. It is necessary then to consider non-thermal 
, mechanisms for desorption. But, if such mechanisms are considered for the production of methanol, they must be considered for all surface 
species. 

Methods. Our gas-grain network of reactions has been altered by the inclusion of a non-thermal desorption mechanism in which the exother- 
micity of surface addition reactions is utilized to break the bond between the product species and the surface. Our estimated rate for this process 
Q ' derives from a simple version of classical unimolecular rate theory with a variable parameter only loosely contrained by theoretical work. 

| Results. Our results show that the chemistry of dark clouds is altered slightly at times up to 10 6 yr, mainly by the enhancement in the gas-phase 
CZ) . abundances of hydrogen-rich species such as methanol that are formed on grain surfaces. At later times, however, there is a rather strong change. 

Instead of the continuing accretion of most gas-phase species onto dust particles, a steady-state is reached for both gas-phase and grain-surface 
^> species, with significant abundances for the former. Nevertheless, most of the carbon is contained in an undetermined assortment of heavy 
• ^ ■ surface hydrocarbons. 

<S Conclusions. The desorption mechanism discussed here will be better constrained by observational data on pre-stellar cores, where a significant 
^ , accretion of species such as CO has already occurred. 
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1. Introduction 

Methanol (CH3OH) is a molecule commonly detected over a 
wide range of conditions in interstellar clouds. In quiescent 
dark cloud regions, it is present in the gas phase with a typi- 
cal abundance of ~1.5 x 10~ 9 «h dSmifh et al.lr2004l) . Gas-phase 
chemical kinetic models have long been successful at reproduc- 
ing this abundance, invoking the radiative association of CH3 
and H2O to form protonated methanol, followed by recombi- 
nation with electrons to produce methanol and atomic hydro- 
gen. The rate coefficient adopted in both the OS U chemical net- 
work (osu.2003) and the UMIST rate99 ratefile dLe Teuff et al 
2000) for the radiative association reaction was estimated to 



be k RA = 5.50 x 10- 12 (77300 r L7 cm 3 s" 1 (for T = 10 - 50 



K) (!Bateslll983l:lHerbstlll985l) . However, the rate c onstant has 
more recently been experimentally determined by iLuca et al 
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(2002); they obtain an upper limit of 2.0 x 10~ 12 cm 3 s _1 for a 
temperature of 50 K, with an uncertainty of 30 K. On the as- 
sumption that this rate is applicable to temperatures of 10 K, 
the canonical gas temperature of quiescent dark cloud regions, 
the measured rate falls short of the estimated value by approx- 
imately 3 orders of magnitude. As outlined by Garrod et al. 
(2006a, hereafter GPCH) such rates are incapable of reproduc- 
ing observed dark cloud methanol abundances. 

In addition, the rate coefficient and branching ratios of the 
dissociative recomb ination of CH^QFTj a nd CD3OD2 were re- 
cently measured by lGeppert et alJ ( 20061) . Whilst their experi- 
ment suggested a rate constant at 10 K a few times larger than 
that previously adopted in the reaction networks, they mea- 
sured a branching fraction for CH3OH of 3 + 2% (6 + 2% for 
CD3OD), as compared to an assumed 50 - 100 % in chemical 
models. The strongest channels produce multiple fragments, 
mainly consisting of separate carbon- and oxygen-groups. The 
combination of these new rates and branching ratios leaves the 
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Table 1. Models 



Table 2. Initial abundances of H 2 and elements (IGraedel et al 
1982). 



Model 


a 


MO 





Ml 


0.01 


M2 


0.03 


M3 


0.1 



gas-phase route incapable of producing the observed methanol 
abundances; peak model values are at least five orders of mag- 
nitude short of the target. Even with an alternative to the radia- 
tive association - a reaction between formaldehyde and proto- 
nated formaldehyde to produce protonated methanol and car- 
bon monoxide - the modelled values are entirely inadequate. 

Hence, in the absence of some heretofore unexplored gas 
phase reaction, we must conclude that grain-surface formation 
of methanol is responsible for gas phase abundances. Surface 
methanol is typically assumed to be formed via the repetitive 
hydrogenation of accreted CO by atomic hydrogen, which is 
comparatively mobile on dust-grain surface s, even at 10 K (see 
Katz et all 1 1999b. Exper i menta l work by IWatanabe & Kouchi 
d2002h and lHidaka et a l. (2004) has demonstrated that this pro- 
cess is efficient on ice surfaces at such low temperatures, in 
spite of activation energies on the order of 1000 K for the re- 
actions H + CO — > HCO and H + H 2 CO -> CH 3 0/CH 2 OH. 
Methanol is indeed detected in icy dust-grai n mantles in in- 



frared absorption towards background stars. iNummelin et al 
(2001) suggest an abundance of < 4 x 10~ 6 n H along the 



line of sight to the b ackground star Elias 1 6. The OSU gas- 
grain chemical code (IHasegawa et al.lll992l) . which uses rate 
equations to model the coupled gas-phase and grain-surface 
chemistry, is capable of producing surface-bound methanol in 
such quantities, hence we would require only around 0.1% of 
total methanol formed in this way to be present in the gas- 
phase to reproduce observed gas-phase abundances. But, cru- 
cially, thermal desorption at 10 K is negligible for species more 
massive than molecular hydrogen; we therefore require some 
other means of returning surface-produced methanol to the gas 
phase. 

Certain non-thermal desorption mechanisms have been im- 
plemented in chemical models; however, only cosmic ray- 
induced heati ng desorption is frequently employed in dark 
cloud models (IHasegawa & Herbstlll993l) . although even this 
mechanism appears too slow to produce appreciable methanol 
desorption. Some models utili s e a d i rect photo-desorption 



mechanism fl_D raine & Salpeterl 11979c lHartquist & Williams 
1990c IWillacv & Millarlll998l) . but the yield per photon, and 
therefore the overall rate, are uncertain. 

Som e work h a s been conducted on other possible mech- 
anisms. IWilliamsl ( 1 19681) first suggested that the stabilisation 
of the excited product of an exothermic grain surface reac- 
tion could result in its deso rption from the sur f ace. Th is idea 
was considered further by Watson & Salpeter ( 1972a, b), but 
no firm desorption fraction was suggested; however, their es- 



Species i 



rii/n H ' 



H 2 


0.5 


riC 


n 1 a 


c + 


7.3(-5) 
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2.14(-5) 
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1.76(-4) 


s + 


8.0(-8) 


Na + 


2.0(-9) 


Mg + 


7.0(-9) 


Si + 


8.0(-9) 


P + 


3.0(-9) 


Cl + 


4.0(-9) 


Fe + 


3.0(-9) 



+ a(b) = ax 10* 



timates were more effic ient than we require in this model. 
Jones & Williams! ( 1984 ) calculated that the products of the re- 
actions O + H — > OH and OH + H — > H 2 should be retained 
on the grains with a probability of between 70 and 100% for 
grain-mantle formation to be efficient in dark clouds. 

Duley & Williams! dl993h suggested that the highly 



exothermic reaction between hydrogen atoms on grain sur- 
faces, forming ro-vibrationally excited H 2 , may impart some 
energy to the grain surface. This would allow the localised 
evaporation of weakly bound molecules like CO in regions 
where H 2 formation had not yet reduced gas-phase atomic hy- 
drogen abundances to their low steady-state levels. 

In this study, we further investigate the non-thermal desorp- 
tion mechanism presented by GPCH, in which the energy of 
formation released by the association reaction of grain-surface 
radicals may break the surface-molecule bond of the product. 
The probability for this to occur is small (on the order of 1 %), 
and therefore most of the product r emains on the grain surface . 
The molecular dynamics study of iKroes & Anderssonl J2006) 
also suggests that this mechanism may be important for pho- 
todissociated water ice. Here we apply the mechanism not only 
to methanol production but to all surface reactions which result 
in a single product. A statistical theory is used to determine the 
probability of desorption, which is dependent on the product 
species and the exothermicity of the reaction. 

We use the OSU gas-grain chemical code to investigate the 
detailed chemistry which takes place both in the gas-phase and 
on the grain-surfaces as a result of the new mechanism. In ad- 
dition, by comparison with observed abundances, we weakly 
constrain the generalised free-parameter, a, which strongly de- 
termines the desorption probability for all qualifying reactions. 

In Section 2 we describe the chemical model and the new 
desorption mechanism. We present the results in Section 3, and 
compare with observations of L134N and TMC-1CP in Section 
4. We present our conclusions in section 5. 
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Fig. 1. Gas phase abundances for models (a) without the new desorption mechanism, MO, and (b) with the desorption mechanism 
at its strongest, M3. 



2. Chemical model 

2.1. The gas-grain code 

To model the chemistry of a dark cloud, we model the gas- 
phase and grain-surface chemistry in tandem using rate equa- 
tions. We allow accretion onto the grain surfaces from the 
gas phase. Surface-bound species may evaporate thermally, ac- 
cording to a Boltzmann law, or by the sporadic heating o f 
the grains by cosmic-ray impacts dHasegawa & Herbstll 1993b . 
which is averaged over time. Surface species may be pho- 
todis sociated, by both the c osmic ray-induced radiation field 
of Prasad & Tarafdarl ( 11983b and by the (heavily extinguished) 
external radiation field. The surface photodissociation rates 
are identical to their g as-phase analogues, as explained by 
Ruffle &Herbs3d2001al) . 

Surface species may react with each other; rates are depen- 
dent on the concentrations, and t he sum of the dif f usion rates 
of the reactants. As explained by lHasegawa et al. ( 1992 ). the 
diffusion rate is defined as the frequency at which the species 
may thermally hop over the barrier between sites, divided by 
the total number of sites on the grain surface. 

2.2. Desorption via exothermic surface reactions 

In addition to the two basic desorption mechanisms, we include 
that introduced by GPCH. For each surface reaction that leads 
to a single product, we apply a branching ratio that this product 
may be ejected into the gas phase, on the assumption that there 
exists a probability that the energy released in the formation 
of the molecule acts to break the surface-molecule bond. We 
assume, in the case of two-product reactions, that the energy is 
lost to lateral translation along the grain surface; hence, in this 
case, no desorption occurs. 

To quantify the probability of desorption, we apply Rice 



Ramsperger-Kessel (RRK) theory (see e.g. iHolbrook et al 



molecular vibrational mode, RRK gives the probability, P, for 
an energy E > Ed to be present in the bond, from a total energy 
F 



1 



(1) 



where Ed is the desorption energy (binding energy) of the 
product molecule, E reac is the energy of formation released in 
the reaction, and s is the number of vibrational modes in the 
molecule/surface-bond system. This number is s — 2 for di- 
atomic species; for all others, s = 3N - 5, where N is the 
number of atoms in the molecule, which is assumed to be non- 
linear. In the hypothetical case where the molecule has no other 
means of energy loss (i.e. to the grain surface), we would ex- 
pect desorption at a rate vP, where v is the frequency of the 
surface-molecule bond (comparable to a Van der Waals bond). 

In practice, we expect energy loss to the surface to be fast, 
so in order to obtain the branching fraction for desorption, we 
must model the competition between the two processes. We 
define a rate at which the total energy is lost to the surface, v s , 
hence the fraction of reactions resulting in desorption is: 



vP 



aP 



v s + vP 1 + aP 



(2) 



where a = v/v s , the ratio of the surface-molecule bond- 
frequency to the frequency at which energy is lost to the grain 
surface. (Whilst the use of a value a + 1 strictly constitutes an 
empirical modification to the pure RRK treatment, such mod- 
ifications are frequently employed to ensure agr eement with 



experi ment in other applications of the theory; see All ain et al 



1996). Modelling the surface-molecule bond as an additional 



(1996), and references therein, for a discussion of this issue). 
In this study we use a generic value of a for all product species. 
Since Ed is normally much less than E reac , P is approximately 
unity and f ~ a, for small a. 

The new mechanism is incorporated into the code such that 
a fraction, /, of the product species in qualifying reactions is 



4 



Garrod et al.: Desorption via exothermic surface reactions 



10 " 


. (a) Gas phase CH 3 0H 


- 


10" 8 






10 -,o 




\Vm2 


10" 12 




MO 


10" 14 


^^^f . . r 


A ; 



Time (yrs) 



Fig. 2. CH 3 OH abundances. 
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Fig. 3. CO abundances. 



desorbed, whilst the rest, (1 - /), remains as a surface-bound 
product. 

G PCH assume a = 0.1 , but this value has b een la belled 
high dPillinjd |2006h . Also, iKroes & Anderssonl d2006l) have 
conducted molecular dynamics simulations of the irradiation 
of water ice with UV photons. They measure the occurrence of 
a number of possible outcomes, including that in which the H 
and OH resulting from photodissociation of H2O recombine. 
From their data, we estimate that ~0.9 % of recombinations re- 
sult in desorption. Using our value of EoiHiO) = 5700 K and 
E reac = 5.91 x 10 4 K, this implies a = 0.012. In order to test 
the effects of the new mechanism, and constrain the value of a, 
we investigate models (see Tabled with various values of a: 0, 
0.01,0.03 and 0.1. 

2.3. Rates and initial conditions 

For this study, we adopt our latest gas-phase chem- 
ical network, osu.2005. The full gas phase rate- 
file and documentation of updates are available at 
http://www.physics. ohio-state. edu/~eric/research.html We 
have updated the grain surface photodissociation rates in 
line with the gas-phase values included in the new network. 
To calculate the various desorption rates, and the diffusive 
reaction rates, we use the binding energies, Ep(i), a nd dif - 
fusion barriers, E B {i), adopted by iGarrod & Herbst (2006), 



corresponding to an amorphous water ice surface. These 
values are typically a little larger than the bare-grain values 
used by GPCH and earlier models (Cuppen & Herbst, in 



prep.) . The review of experimental evidence by iKatz et al 



1999h indicates that quantum tunnelling through diffusion 



barriers is inefficient, even for atomic hydrogen, therefore all 
diffusive rates are based o n a thermal hopping rate. We treat the 



so-called modified rates (Caselli et al. 1998; Stantchevaet al 



2001) 



including reactions with ac t ivation energy barriers, in 
the same way as IGarrod & Herbst (120061) : all surface atomic 
hydrogen reaction rates may be modified. At 10 K, surface 
reactions are domina ted by H-addition ; other species are much 
less mobile. As per IGarrod & H erbst (2006), we assume an 
activa tion energy, Ea = 2500 K d Woonll2002l: iRuffle & Herbst 
l200d) . for both of the reactions H + CO -> HCO and H + 
H 2 CO -> CH3O. 



Initial abundances correspond to the so-called low metal 



abundances of lGraedel et al.1 (119821) : see table [2] Physical con- 
ditions remain constant throughout the model, with a density 
n H — 2 x 10 4 cur 3 , gas and grain temperatures T K — T„ — 10 
K, and a visual extinction Ay = 10. The cosmic-ray ionisation 
rate is set to the canonical value, £ = 1.3 x 10~ 17 s~'. 
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Fig. 4. Fraction of carbon contained by most important carbon- 
bearing species for model M0. Surface-bound species are rep- 
resented by (s). 

3. Results 

Figures QJ & b show the calculated abundances of a selection 
of gas phase species for the two extremes of the model: without 
the new desorption mechanism, model M0, and with the new 
desorption mechanism at its strongest, i.e. a = 0.1, model M3. 
These plots show the results of using our amorphous water ice 
surface, rather than the silicaceous/carbonaceous surface used 
in previous dark cloud gas-grain models (GPCH). Most species 
are more tightly bound to the surface, resulting in smaller late- 
time abundances. Indeed, in model M0, where the strongest 
desorption mechanism is that caused by cosmic ray-induced 
heating, desorption is now insignificant for all but the lightest 
species. For example, the new CR-induced desorption rate for 
methanol is negligible, using the measured binding energy of 
5530 K dCollings et al.ll2004 . 

The effects of the new mechanism at late times are quite 
stark - model M3 does not display the total gas phase deple- 
tion of heavy species seen in previous models, but reaches an 
approximate steady state by a time t — 10 8 yr. Models Ml 
(a = 0.01) and M2 (a = 0.03) also reach steady state, but abun- 
dances are proportionately lower, by approximately lOx and 
3x, respectively. Thus, for the majority of gas phase species, 
the strength of the effect of the new mechanism is quite pre- 
dictable at very late times. 

At times before ~1 Myr, the majority of species do not 
exhibit very significant variations between models, although 
many species do show a somewhat greater peak value (at 
around 10 6 - 10 7 yr) for a + 0. Since, at the grain temperature 
of T g = 10K, hydrogenation is the dominant chemical reaction, 
the effect of the new desorption mechanism is most strongly 
represented in the gas phase by hydrogen-bearing species. Most 
strongly affected are the multiply-hydrogenated species, since 
their partially hydrogenated precursors are also ejected into the 
gas phase, where their hydrogenation may also be completed. 
The peak H2O abundance is raised due to the formation and 
partial desorption of large amounts of OH and water from the 
grain surfaces; however, water is still the dominant ice compo- 
nent, and the build-up of water ice mantles is not significantly 
affected at any time by the new mechanism. CO is also slightly 



enhanced at its peak, partially due to the desorption of OH, 
which facilitates reaction with ionised carbon. Naturally, the 
desorption of OH tends to increase the abundance of O2 (not 
shown), via reaction with atomic oxygen. H2S is very strongly 
enhanced, and is a special case, as discussed in section 3.2. 

3.1. Selected carbon- and oxygen-bearing species 

The peak methanol abundance produced in GPCH with param- 
eter a = 0.1 is somewhat higher than observations suggest, typ- 
ically a few x10~ 9 «h- In model M3, this effect is mitigated to 
some degree by the new stronger binding energies; the resultant 
peak abundance, 1.5 x 10~ 8 /ih, is high, but just within an order 
of magnitude of the typical levels. Figure [2] shows gas-phase 
and grain-surface methanol abundances for each model. With 
a + 0, gas-phase methanol abundances scale approximately 
with a; model Ml achieves 1.8 x l0^ 9 n H at peak. The range of 
peak values for models Ml - M3 confirms the success of the 
new desorption mechanism in producing appropriate levels of 
gas-phase CH3OH. 

The surface methanol abundance is largely unaffected by 
the new mechanism until very late times. Interestingly, at these 
late times the surface CH3OH abundance increases with in- 
creasing a, so that a larger desorption parameter does not result 
in lower surface methanol levels; indeed they are strongly en- 
hanced. As can be seen from Fig. [3] the behaviour of CO, the 
other main repository of grain-surface carbon, is very similar 
to that of methanol. Figures [2] and [3] also show commensurate 
increases in gas phase CO and CH3OH. What is happening at 
late times? 

To answer this question, it is first necessary to consider the 
major carbon-bearing species as a function of time. Figure [4] 
shows these dominant species, for model M0 (a = 0). We see 
that in this model (and in models Ml - M3), most of the carbon 
does not end up in CO or CH3OH but at late times goes through 
a period in which solid methane dominates before ending up 
in the form of an assortment of larger surface-bound hydro- 
carbon molecules. The state of hydrogenation of such surface- 
bound hydrocarbons is poorly constrained in our model. The 
long timescales required (? >~10 7 yr) may of course make the 
attainment of large hydrocarbon abundances impossible, due to 
dynamical considerations. 

The conversion of surface CO and CH3OH into methane 
via a variety of processes is depicted in Fig. [5] The methane 
itself is largely formed by hydrogenation of the methyl 
(CH3) surface radical. Without the new desorptive mecha- 
nism, cosmic -ray photodissociation of CH3OH is the primary 
formation route for CH3 at later times; the methanol CR- 
photodissociation branches that produce CO and H2CO quickly 
result in re-hydrogenation, producing no net effect. Cosmic-ray 
photodissociation of methane also quickly leads to its own re- 
formation. 

With the new mechanism activated, the conversion of CO 
and CH3OH to CH4 occurs during continual recycling between 
gas and grain surface - the direct dissociation of methanol 
into CH3 is no longer dominant. These routes rely heavily on 
the initial dissociation of surface methanol by cosmic-ray in- 
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duced photons and desorption via the new mechanism, leading 
to a variety of gas-phase species that undergo an assortment 
of ion-molecule and neutral-neutral processes or re-accrete 
onto grains. The gas-phase processes eventually produce CH3 
among other species. This radical then can accrete onto the sur- 
face as a precursor for surface methane and, via the new desorp- 
tive mechanism, gas-phase methane. Alternatively, it can react 
with atomic C to form acetylene (C2H2) and H atoms in the 
gas, leading to more complex hydrocarbons. In the gas phase, 
methane can also be converted into larger hydrocarbons by a 
variety of processes. 

Whilst without the new mechanism, most surface methanol 
is ultimately channelled into CH4 at late times, the new mech- 
anism allows carbon hydrides to return to the gas phase, where 
they may be converted into CO or its hydrides, maintain- 
ing modest levels of CO, formaldehyde and methanol in both 
the gas phase and on grain surfaces. The continuous CR- 
photodissociation of surface CH4 (and other hydrocarbons) 
keeps this cycle active. These processes are also shown in 
Fig.S 

We may comp are the late-time beh aviour of our new model 
with the results of Willacv et al. (1994). They investigated the 
mechanism of IDuley & Williams! (119931) . whereby the forma- 
tion of H2 releases energy into the grain surface, allowing lo- 
calised evaporation of weakly bound molecules, most notably 
CO. Whilst we use a much larger gas-phase and grain-surface 
reaction network, whose rates have also evolved somewhat in 
the intervening period, we may comment broadly on the most 
obvious differences. Whereas we find that complete freeze-out 
onto grains is prevented by our new mechanism, theirs has the 



effect of merely retarding freeze-out. This is due to our in- 
clusion of grain-surface photo-destruction processes, which al- 
lows re-hydrogenation of the products and the associated evap- 
oration via the new mechanism. The mechanism of Willacy et 
al. allows significant re-injection of mantles into the gas phase 
only while the accreted species are small, with low binding en- 
ergies, such that the formation energy of reaction H + H — > H2 
may desorb them. When CO is converted in the gas phase into 
OH and H2O, and CH4 and C2H2, and subsequently accretes, 
the mechanism is no longer effective and freeze-out takes over. 
Note also that surface CO is left unprocessed in their model, 
which tends to prolong the efficiency of their mechanism in re- 
turning carbon and oxygen to the gas phase. 

3.2. Sulphur-bearing species 

The behaviour of sulphur is a long-running problem in inter- 
stellar chemistry. For example, gas-phase SO and CS are typi- 
cally over-produced. Previous gas-grain models have predicted 
large quantities of H2S on grain surfaces; however, no solid 
phase H2S has been detected to date in the interstellar medium. 
As in this model, the initial elemental sulphur abundance is 
therefore usually reduced by around 2 orders of magnitude on 
the diffuse cloud abundance. 

Figures|6h & b show the gas-phase and grain-surface abun- 
dances of H2S for each model. Model MO produces too much 
surface H2S (i.e. a significant proportion with respect to H2O 
ice), and too little in the gas phase. However, the introduction 
of the new mechanism drastically reduces the abundance on the 
surface, and greatly increases the gas-phase level, now much 
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Fig. 6. H 2 S abundances. 



more in li ne with typical obs erved values (e.g. 8x10 I0 «(H2) 
in L134N. IOhishi et aljl992h . from around 3 x 10 4 yr onwards. 
Surface abundances are around 10~ 5 of the H2O ice abundance, 
even for model Ml, acceptably low to agree with observations. 
The inclusion of the new mechanism, for the range a = 0.01 - 
0.1, also produces abundances of gas-phase CS and SO in good 
agreement with observational values at appropriate times, al- 
beit using a much depleted initial abundance of atomic sulphur. 

The unusual strength of the effect on H2S is due to the 
hydrogen abstraction reaction H + H2S — > HS + H2. H2S 
is formed by repetitive hydrogenation of atomic sulphur on 
the grain surfaces. Because the abstraction reaction involves 
atomic hydrogen, and has only a fairly small activation energy 
(860 K), it is relatively fast. In model M0, the ratio of HS to 
H2S achieves a quasi-equilibrium, according to the relative hy- 
drogenation and abstraction rates, resulting in the steady for- 
mation of H2S in tandem with H2O. For models Ml - M3, 
although the abstraction reaction itself does not (in our model) 
directly lead to desorption via the new mechanism, the resul- 
tant HS is quickly hydrogenated and this reaction does lead to 
desorption. Sulphur, in its atomic and hydrogenated forms, is 
stuck in this fast loop, and so the new mechanism is very effi- 
cient at "syphoning off" sulphur back into the gas phase. The 
result is that at times from around 10 4 - 3 x 10 6 yr, gas phase 
atomic sulphur becomes the dominant sulphur-bearing species 
as the desorbed H2S is broken down. Through the rapid H addi- 
tion/abstraction rates, the process of removal is efficient enough 
to ensure a minimal a-dependence for gas-phase H2S abun- 
dances; no significant build-up of grain-surface sulphur occurs 
for models Ml - M3. In all models, surface-bound H2CS be- 
comes the dominant form of sulphur, at times later than ~5 
Myr. 

The process of H2S formation and re-formation on grains 
may have implications for the degree of deuteration detected 
in this species. The speed of the hydrogenation and abstraction 
reactions would allow a quasi-equilibrium between deuterated 
and undeuterated forms to arise quickly; deuterium fractiona- 
tion should therefore be highly sensitive to the branching ra- 
tios of the abstraction reactions (e.g. for H + HDS, the ratio of 
products H2 + DS, versus HD + HS), as well as the specific 
activation energy barriers for each reaction. 



10" 
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M0 - 
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Fig. 7. NH3 abundances. 



3.3. Nitrogen-bearing species 

The effects of the new mechanism on nitrogen surface chem- 
istry are not very substantial - the majority of nitrogen still ends 
up as NH3 o n the grains, a result in some disagreement with 
observation (IWhitteM I2003I) . When ammonia is broken down 
on grain surfaces by cosmic ray-induced photodissociation into 
NH2 or NH, these products can be hydrogenated, to ammonia 
and NH2 respectively, and can be desorbed according to the 
new mechanism. However, the gas-phase chemistry also tends 
to favour the formation of ammonia from NH2, and the ammo- 
nia may ultimately re-accrete. Even for the a = 0.1 case, the 
primary form of surface nitrogen is unchanged. 

An important effect is the raised abundance of gas-phase 
NH3 at early times (see Fig. 13, due to the very moderate frac- 
tions released from the grain surfaces. At nearer to 1 Myr, the 
difference among models tails off, and the peak value, at t ^ 5 
Myr, is unaffected. Clearly at times of 10 7 yr or greater, there 
is a very large disparity between M0 and the models which 
include the new mechanism. Hence, the new desorption mech- 
anism has the effect of making NH3 more of an "early-time" 
species, although the effect is very significant only for the 
largest a-value. 

The cyanopolyynes included in the code (HC3N, HC5N, 
HC7N & HC9N) are all enhanced at both their early-time peak 
and late-time peak; the enhancement is an order of magnitude 
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for model M3. This effect arises primarily from the injection of 
hydrocarbons into the gas phase, allowing increased reaction 
with atomic nitrogen. 

4. Comparison with observations in quiescent 
dense cores 

The analysis of GPCH included a comparison of the model 
with the observational results for TMC-1C P. This was done 
in the same fashion as Smit h et al. (2004), simply by com- 
paring the number of species with computed abundances that 
fell within one order of magnitude of the detected level. This 
method is helpful in comparing models which are very differ- 
ent, i.e. gas phase vs. gas-grain, however it is less useful for 
the comparison of the range of models which we exhibit here, 
which are mostly quite similar. It is also harder to distinguish 
the time of best fit for a particular model. The problem arises 
from setting such a strict criterion for a "success"; inevitably, 
some species fall close to the order of magnitude limit, and so 
similar models can produce wildly different levels of "success". 
We therefore require a quantitative means of evaluation of the 
success of the model in repr oducing observa t ions. 

As a different approach, Wa kelam et al. compared 
observations with the results of a gas-phase model taking into 
account both observational and theoretical error bars, the latter 
computed using statistical methods based on uncertainties in 
rate coefficients. Agreement for each species was defined by 
an overlap between error bars, and a logarithmic distance of 
disagreement was computed for each species in the absence 
of overlap. Due to the large number of parameters in the gas- 
grain code, a full analysis of the error propagation would be 
very difficult to achieve, and has not so far been attempted. 

Here we present a new means of comparison, whereby we 
assign a level of confidence in the agreement of each computed 
abundance (at a particular time) with the observed value, hav- 
ing designated it with a generic observational error. We con- 
struct a log-normal distribution about each observational value, 
and identify its defining standard deviation, cr, with an appro- 
priate error factor on the observed value. In a similar way as 
one might determine the confidence level of a spectroscopic 
detection being distinct from the mean baseline, we conversely 



define our "confidence", k,, that the calculated value, X), for 
species i is associated with the observed value, X b S j, thus: 



Ki — erfc 



log (Xd- log (X, 



obs 



y/2o- 



(3) 



where erfc is the complementary error function {erfc — 1 - erf); 
Ki ranges between zero and unity. The term "confidence" is con- 
tained within quotation marks to signify that we are not dealing 
with a rigorous statistical analysis. For our analysis, we define 
cr = 1, hence 1 standard deviation corresponds to one order of 
magnitude higher or lower than the observed value. Therefore, 
a calculated value which lies 1 order of magnitude from X i c ,i, s 
has a confidence level a-, = 0.317, whilst a value 2 orders of 
magnitude from X b s j has a confidence level = 0.046. For 
each time in the model, we take a mean of the individual confi- 
dence levels of every species, allowing us to define the overall 
confidence in a match with observations for each model at each 
time. 

Other means of comparing quantitatively with observations 
were considered; for example, taking the sum or mean average 
of | log Xi — log X bs.i\ ov e r all s pecies, in a similar approach 



to that of Wakelam et al. (2006), to provide a parameter (the 
mean "deviation") which becomes smaller with a closer match. 
However, this method tends to skew results to some degree 
when there are "large" deviations from the observed values. 
For example, if in one model a particular species is 3 orders 
of magnitude from the observed value, and in another model, 
that species is 6 orders of magnitude away, the influence of 
this species on each mean will be large, even though in both 
cases the actual value is so far away as to be no longer cred- 
ible. In such a case, we would suggest that the model itself is 
flawed (as regards that species) and that no amount of "tun- 
ing" can reasonably be expected to produce the right result (for 
the right reasons). But this would not necessarily make either 
model catastrophically wrong, in spite of what the comparison 
parameter might su ggest. The "yes or no" method employed by 
Smith et al.l ([2004) naturally discards such outliers. Our pre- 
ferred method, as explained above, allows us to define, in a 
quasi-statistical manner, the extent to which variations in cal- 
culated values are important when the values themselves are 
far from the target. Whilst our method is no more statistically 
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rigorous than the alternatives adumbrated above, it provides a 
convenient and self-consistent means of analysis that is easily 
automated, and which allows quantitative comparison between 
different models. 

We now compare the results with detected abundances ob- 
tained for L134N and TMC-1CP (the so-called cyanopolyyne 
peak). Given in Tables [3] and |H these values are taken from a 
number of different stu dies, but corresp ond to those collate d by 
Wakelam et all d2006l) . for L134N, and lSmith et ail d2004 . for 
TMC-1CP. To those species that have only an upper or lower 
limit, we assign a confidence level of unity if the calculated 
value is lower than the upper limit, or greater than the lower 
limit, respectively; they are otherwise treated in the same way. 



4.1. L134N 

Figure [8] shows the level of confidence in the match between 
each model and the 41 observational values for L134N, for 
times 10 4 - 10 s yr. Clearly the models Ml - M3 are an im- 
provement over the values obtained without the new mecha- 
nism. The best match at any time is obtained for model M2 
(a = 0.03), at a time of ~3 Myr; however, each of models 
Ml - M3 reaches its peak at around this time, and the differ- 
ences in confidence are very small. The model M2 peak is also 
the widest, spanning around 1 Myr with a ~69% confidence 
level. Each of the models shows a confidence of ~65 % over a 
wider time interval of t — 10 5 - 3 x 10 6 yr. This corresponds 
to less than a factor of 3 deviation, on average, between mod- 
elled and observed abundances. Model M3 in particular has an- 
other strong peak, at ~10 5 yr, but M2 also shows a local peak 
at around this time. An age for L134N of anywhere between 
0. 1 and 3 Myr may therefore be plausible according to mod- 
els Ml - M3; methanol is still adequately produced at all such 
times. Whilst M3 also has a very high agreement at 10 8 yr, on 
dynamical grounds we should consider a match at this time to 
be highly dubious. 

Clearly, whilst the introduction of the new mechanism 
much improves the agreement with observed abundances ac- 
cording to our method of analysis, variation of the parameter a 
in the 0.01-0. 1 range does not strongly alter it, although higher 
values make an earlier match more of a possibility. The mean 
confidence level appears therefore to give us only limited scope 
for the constraint of the a-parameter. Comparing the mean con- 
fidence level approach with a method that utilizes the mean of 
the log of the deviation, we find the latter to yield similar times 
of best agreement, but to suggest that M3 (as opposed to M2) 
achieves the best match with observations. The difference is 
marginal, though, and serves to demonstrate the difficulty in 
constraining the value of a. 

Table [3] shows the abundances of molecules observed in 
L134N and the calculated values for each model at its own 
peak at ~2 Myr. Values which vary from the detected value 
by more than one order of magnitude are printed in boldface. 
Values which exceed the detected level by this much are also 
set in italics. 

Aside from SO, H 2 S and O2, all of the highlighted dis- 
crepancies occur in the abundances of carbon-bearing species. 



Indeed, most of these occur for species containing more than 1 
carbon atom. The larger cyanopolyynes are not well matched at 
the peak times, typically falling quite short of observed values. 
At later times (close to 10 Myr), however, these abundances 
are fairly well matched with observations. CN is badly over- 
produced in all models. Hydrocarbon abundances are typically 
raised for higher a-values, although, except for C3H, they are 
usually still a fair match. The best-matched time for larger a 
also tends to under-represent atomic C, though not by more 
than 1 order of magnitude. For the most part, those species 
with observed abundances not well-reproduced by model M0 
are either much improved upon with all of models Ml - M3, 
or barely improved at all - much as Fig.[8]represents. Very few 
species are actually worsened in their agreement by the new 
mechanism. The best match to observations with the order-of- 
magnitude criterion is still model M2, with an age of around 2 

- 3 Myr, the same as determined by the mean confidence level. 
However, earlier matches, to as short an age as 10 5 yr, are plau- 
sible, and methanol is still adequately produced at such times. 

It is worth comparing the success of the gas-grain code 
with that of our most recent gas-phase code. lSmith et al. (2006) 
achieved 73% of calculated abundances falling within 1 order 
of magnitude of the observed value, for L134N, at time 10 5 
yr. Our best model achieves 88% (36/41 species), at ~2 Myr. 
Clearly, our new gas-grain models improve agreement with ob- 
servations over a purely gas-phase analysis, and this is true over 
the entire range of ~10 5 - few xlO 6 yr. However, in the case 
where a = 0, model M0, a comparable level of agreement with 
the gas-phase treatment (70%) is achieved. 

4.2. TMC-1CP 

Figure [9] shows the mean confidence levels for each model in 
comparison to 52 species detected in TMC-1CP. The low mean 
confidence obtained at all times is apparent, and corresponds 
to a factor of 5 - 10 deviation on average, significantly higher 
than for L134N. Nevertheless, we can clearly see that greater 
a-values produce a better match with the TMC-1 values. Also, 
there are two peaks, but the later-time peak is much stronger, 
especially for greater values of a. This stands in contrast to the 
commonly-held view of TMC-1 as a relatively young object 
dHartquist et al.ll200ll) . 

Table |4] shows the observed and computed best-time abun- 
dances for each molecule in the observational data set. We can 
see immediately that most molecules that fail to be reproduced 
fail to some degree in almost every model. As for L134N, the 
majority of these are carbon-bearing species. Some of these 
species achieve an acceptable match with observations using 
model M3 (a = 0.1), however, other carbon-bearing species 
then attain levels which are too high. For other species, model 
M3 merely improves upon the other models, without getting 
reliably close to the observed value. In the case of models M0 

- M2, we might also argue that the fact that so many values 
fall below the observed values, yet so few lie above, makes 
it unlikely that the discrepancies are purely a statistical effect; 
however, model M3 is well-balanced in this respect. Test runs 
using values yet higher than a = 0.1 produce an ever better 
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Table 3. Observational values of fractional abundances in L134N, and corresponding values at the "best fit" times for each model. 



Species 


N(i)/N(H2) '' 2 




n(0/«(H 2 ) 13 






Observed 


MO 


M 1 


M2 


M3 




(L134N) 


t = 2.6 Myr 


/ = 1.9 Myr 


1 = 2.8 Myr 


( = 3.1 Myr 


C 


>1.0(-6) 4 


6.2(-7) 


8.1(-8) 


5.3(-7) 


2.0(-7) 


CH 


1 .0(-8) 5 


6.8(-9) 


3.2(-9) 


7.7(-9) 


7.7(-9) 


C 2 H 


<5.0(-8) 5 


8.5(-9) 


5.5(-9) 


9.3(-9) 


1.1 (-8) 


C 3 H 


3.0(-10) 5 


6.8(-9) 


1.5(-9) 


8.6(-9) 


6.9(-9) 


C 3 H 2 


2.0(-9) 5 


1.0(-10) 


3.6(-ll) 


2.2(-10) 


3.9(-10) 


C1H4 


<1.2(-9) 5 


4.0(-ll) 


2.4(-10) 


3.5(-10) 


1 -4(-9) 


C 4 H 


1.0(-9) 5 


3.9(-9) 


1.1 (-9) 


4.7(-9) 


5.3(-9) 


CN 


8.2(-10) 5 


3.3(-8) 


1.4(-8) 


4.3(-8) 


4.8(-8) 


HCN 


1 .2(-8) 6 


4.5(-8) 


2.0(-8) 


6.9(-8) 


8.3(-8) 


HNC 


4.7(-8) 6 


4.0(-8) 


1.9(-8) 


6.3(-8) 


7.5(-8) 


HiCN + 


<3.1(-9) 5 


8.3(-10) 


3.6(-10) 


1.1 (-9) 


1.2(-9) 


CHiCN 


<1.0(-9) 5 


3.3(-10) 


1.0(-10) 


6.1(-10) 


8.8(-10) 


CHiCN 


<1.0(-9) 5 


1 . 1 (-10) 


2.3(-ll) 


2.0(-10) 


2.4(-10) 


CiN 


<2.0(-10) 5 


4.3(-10) 


1.0(-10) 


5.5(-10) 


5.5(-10) 


C 3 H 3 N 


<1.0(-10) 5 


1.6(-12) 


7.7(-13) 


5.7(-12) 


2.3(-ll) 


HC 3 N 


8.7(-10) 6 


2.2(-10) 


1-1(10) 


4.3(-10) 


8.5(-10) 


HC,N 


1.0(-10) 5 


2.2(-ll) 


7.5(-12) 


4.6(-ll) 


1.2(-10) 


HC 7 N 


2.0(-ll) 5 


S.9(-13) 


1.4(-13) 


1.6(-12) 


5.1(-12) 


CO 


8.0(-5) 5 


1.2(-5) 


2.2(-5) 


1.4(-5) 


2.3(-5) 


HCO+ 


1.0(-8) 6 


2.1(-9) 


3.0(-9) 


2.3(-9) 


3.0(-9) 


H 2 CO 


2.0(-8) 5 


2.1(-8) 


1.3(-8) 


3.7(-8) 


7.8(-8) 


CH3OH 


3.7(-9) 6 


6.0(-13) 


1.7(-9) 


4.1 (-9) 


1.9(-8) 


HCOOH 


3.0(-10) 5 


1.5(-11) 


15(-10) 


3.0(-ll) 


1.8(-10) 


CH 2 CO 


<7.0(-10) 5 


2.9(-10) 


3.9(-10) 


3.2(-10) 


4.1(-10) 


CH,CHO 


6.0(-10) 5 


1.3(-13) 


5.6(-13) 


1.3(-12) 


7.0(-12) 


C3O 


<5.0(-ll) 5 


l.O(-ll) 


7.2(-12) 


1.4(-11) 


2.1(-11) 


H ; S 


8.0(-10) 5 


1.6(-11) 


2.9(-9) 


3.1(-9) 


3.1 (-9) 


SO 


3.1 (-9) 6 


6.9(-ll) 


3.2(-9) 


1.5(-9) 


4.4(-9) 


SO ; 


<1.6(-9) 6 


8.9(-13) 


7.6(-ll) 


3.2(-ll) 


3.3(-10) 


CS 


1.7(-9) 6 


6.8C-10) 


6.1 (-9) 


1.3(-8) 


1.9(-8) 


HCS+ 


6.0(-ll) 5 


2.7(-12) 


2.1(-11) 


5.0(-ll) 


6.3(-ll) 


H 2 CS 


6.0(-10) 5 


2.0(-10) 


1.1 (-9) 


5.K-9) 


8.4(-9) 


C 2 S 


6.0(-10) 5 


2.5(-ll) 


l.S(-10) 


3.5(-10) 


2.6(-10) 


C 3 S 


<2.0(-10) 5 


3.7(-12) 


1.4(-11) 


5.6(-ll) 


4.6(-ll) 


ocs 


2.0(-9) 5 


4.0(-12) 


6.2(-ll) 


7.5(-ll) 


1.4(-10) 


NH 3 


9.1(-8) 5 


1.7(-7) 


1.2(-7) 


2.3(-7) 


2.7(-7) 


N 2 H + 


6.8(-10) 6 


2.3(-9) 


2.1 (-9) 


2.5(-9) 


2.2(-9) 


NO 


6.0(-8) 5 


4.4(-8) 


9.1(-8) 


5.9(-8) 


1.4(-7) 


OH 


7.5(-8) 5 


3.6(-8) 


4.9(-8) 


4.9(-8) 


1.1 (-7) 


H 2 


<3.0(-7) 7 


4.5(-8) 


8.6(-8) 


9.K-8) 


2.6(-7) 


O2 


<1.7(-7) 8 


1.8(-7) 


2.6(-6) 


2.0(-7) 


8.8(-7) 


Matches 9 /41 




29 


33 


36 


35 



1 a{b) - ax 10 fo ; 2 Observed values collated by Wakelam et al. (2006); 3 Boldface indicates a 
theoretical value different by more than 1 order of magnitude from the observed value. Plain 
boldface indicates too low a value; italic boldface indicates too high a value; 4 Stark et al. (1996); 
5 Ohishietal. (1992) ; 6 Dickens et al. (2000) ; 7 Snell et al. (2000) ; 8 Pagani et al. (2003) 
9 Agreement with observational value, to within 1 order of magnitude. 



match, however, such values would be very difficult to justify 
on physical grounds, and the fact that their agreement peaks at 
even later times makes the match less plausible. The implica- 
tion is not therefore that M3 is a significantly better model of 
TMC-1CP, but that each model has failed, model M3 having 



failed only to a lesser degree. Dismissing the late-time peak on 
chemical and dynamical grounds would therefore suggest an 
optimum age of around 1 - 3 x 10 s yr, in keeping with typical 
estimates. 
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Table 4. Observational values of fractional abundances in TMC-1CP, and corresponding values at the "best fit" times for each 
model. 



Species 


iV(i)/W(H 2 ) 1J 




«(/)/n(H 2 ) 1J 






Obsei"ved 


MO 


Ml 


M2 


M3 




(TMC-1CP) 


/ = 3.2 Myr 


t = 3.1 Myr / 


= 3.3 Myr 


t = 6.0 Myr 


CH 


2(-8) 


8.5(-9) 


8.7(-9) 


9.0(-9) 


8.7(-9) 


c 2 


5(-8) 


4.4(-8) 


3.8(-8) 


5.3(-8) 


2.4(-7) 


C 2 H 


2(-8) 


8. K-9) 


8.2(-9) 


8.4(-9) 


7.3(-9) 


C 3 H 


l(-8) 


1.2(-8) 


1.2(-8) 


1.5(-8) 


3.0(-8) 


C 3 H 2 


l(-8) 


1.9(-10) 


2.3(-10) 


4.0(-10) 


2.6(-9) 


C3H4 


6(-9) 


7.4(-ll) 


1.7(-10) 


4.2(-10) 


6.4(-9) 


C 4 H 


9(-8) 


5.9(-9) 


5.6(-9) 


7.7(-9) 


4.9(-8) 


C 4 H 2 


K-9) 


2.3(-9) 


2.2(-9) 


3.2(-9) 


2.3(-S) 


C 5 H 


6(-10) 


7.5(-10) 


7.7(-10) 


1.1 (-9) 


4.9(-9) 


CH 3 C 4 H 


4(-10) 


4.9(-ll) 


5.5(-ll) 


9.8(-ll) 


1.7(-9) 


C 6 H 


2(-10) 


2.0(-10) 


2.0(-10) 


3.3(-10) 


3.2(-9) 


QH 2 


5(-ll) 


5.4(-ll) 


6.0(-ll) 


l.l(-lO) 


1.7(-9) 


CN 


5(-9) 


5.0(-8) 


5.1(-S) 


6.4(-S) 


1.2(-7) 


HCN 


2(-8) 


7.8(-8) 


8.5(-8) 


1.1 (-7) 


2.7(-7) 


HNC 


2(-8) 


7.1 (-8) 


7.6(-8) 


1.0(-7) 


2.5(-7) 


H 2 CN + 


2(-9) 


1.5(-9) 


1.5(-9) 


1.8(-9) 


3.2(-9) 


CH 2 CN 


5(-9) 


7.0(-10) 


7.9(-10) 


1.1 (-9) 


3.6(-9) 


CH 3 CN 


6(-10) 


2.6(-10) 


2.9(-10) 


4.1(-10) 


1.1 (-9) 


C 3 N 


6(-10) 


7.0(-10) 


7.2(-10) 


9.5(-10) 


3.8(-9) 


HNC 3 


6(-ll) 


3.9(-ll) 


4.K-11) 


5.8(-ll) 


2.8(-10) 


HC 2 NC 


5(-10) 


3.2(-ll) 


3.5(-ll) 


5.3(-ll) 


3.6(-10) 


C 3 H 2 N + 


1(-10) 


8.7(-12) 


8.7(-12) 


1.2(-11) 


5.9(-ll) 


C 3 H 3 N 


4(-9) 


4.8(-12) 


6.1(-12) 


1.4(-11) 


3.0(-10) 


CH 3 C 3 N 


8(1 


8.4(-13) 


9.8(-13) 


1.6(-12) 


6.8(-12) 


HC 3 N 


2(-8) 


4.7(-10) 


5.2(-10) 


8.5(-10) 


6.2(-9) 


HC 5 N 


4(-9) 


4.8(-ll) 


5.1(-11) 


9.8(-ll) 


2.2(-9) 


HC7N 


K-9) 


1.7(-12) 


1.9(-12) 


4.3(-12) 


1.8(-10) 


HC9N 


5(-10) 


1.6(-13) 


1.7(-13) 


4.2(-13) 


4.0(-ll) 


CO 


8(-5) 


7.6(-6) 


8.7(-6) 


1.0(-5) 


9.5(-6) 


HCO+ 


8(-9) 


1.5(-9) 


1.6(-9) 


1.7(-9) 


1.0(-9) 


H 2 CO 


5(-8) 


2.5(-8) 


2.9(-8) 


3.8(-8) 


6.2(-8) 


CH 3 OH 


3(-9) 


1.0(-12) 


1.1 (-9) 


3.7(-9) 


1.1 (-8) 


HCOOH 


<2(-10) 


5.8(-12) 


9.6(-12) 


1.8(-11) 


8.5(-ll) 


CH 2 CO 


6(-10) 


1.8(-10) 


2.1(-10) 


2.3(-10) 


2.3(-10) 


CH 3 CHO 


6(-10) 


1.5(-13) 


4.4(-13) 


1.3(-12) 


1.7(-11) 


C 2 


6(-ll) 


4.6(-12) 


5.5(-12) 


7.7(-12) 


8.9(-12) 


c 3 o 


l(-lO) 


8.7(-12) 


10(11) 


1.6(-11) 


1.3(-10) 


H 2 S 


<5(-10) 


1.2(-11) 


3.0(-9) 


3.0(-9) 


7.8(-10) 


SO 


2(-9) 


3.0(-ll) 


8.1(-10) 


1.2(-9) 


1.4(-9) 


SO: 


<l(-9) 


2.8(-13) 


9.9(-12) 


2.2(-ll) 


1.9(-10) 


cs 


4(-9) 


6.9(-10) 


1.5(-8) 


2.0(-8) 


2.0(-8) 


HCS+ 


4(-10) 


3.0(-12) 


6.1(-11) 


7.6(-ll) 


6.2(-ll) 


H 2 CS 


7(-10) 


2.6(-10) 


5.9(-9) 


8.0(-9) 


8.6(-9) 


c : s 


8(-9) 


2.3(-ll) 


4.6(-10) 


5.3(-10) 


3.9(-10) 


C 3 S 


K-9) 


3.7(-12) 


7.4(-ll) 


9.K-11) 


9.2(-ll) 


ocs 


2(-9) 


2.1(-12) 


5.1(-11) 


6.K-11) 


2.8(-ll) 


NH-, 


2(-8) 


2.3(-7) 


2.4(-7) 


2.7(-7) 


3.6(-7) 


N 2 H + 


4(-10) 


2.3(-9) 


2.4(-9) 


2.4(-9) 


1.4(-9) 


NO 


<3(-8) 


2.8(-8) 


3.4(-8) 


4.K-8) 


9.4(-8) 


OH 


2(-7) 


2.6(-8) 


3.1 (-8) 


4.3(-8) 


1.6(-7) 


H 2 


<7.0(-8) 


3.4(-8) 


5-K-8) 


9.0(-8) 


5.0(-7) 


o 2 


<7.7(-8) 


4.0(-8) 


5.5(-8) 


6.4(-8) 


8.4(-8) 


Matches 4 /52 




30 


34 


36 


37 



1 a{b) = a X 10''; 2 Observed values collated by Smith et al. (2004); 3 Boldface indicates a 
theoretical value different by more than 1 order of magnitude from the observed value. Plain 
boldface indicates too low a value: italic boldface indicates too hisrh a value. 
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Table 5. Observed ice composition towards Elias 16, expressed 
as a percentage of H2O abundance. 1 



Species 


Elias 16 


Model M2 






2 x 10 5 yr 


2.8 x 10 6 


H 2 » 


1.28(-4) 


3.47(-5) 


1.50(-4) 


CO 3 


26 


7.4 


0.3 


co 2 4 


20 


4(-2) 


0.3 


CH4 5 




20.7 


28.5 


CH30H 5 


< 3 


6.4 


10.5 


H 2 co 5 




5.7 


1.0 


H 2 2 5 


< 5 


3(-4) 


5(-5) 


OCS 5 


< 0.2 


9(-13) 


6(-14) 


NH 3 5 


< 9 


12.6 


7.3 



1 a(b) -ax 10 b ; 2 H 2 Q fraction al abundance; 3 re f . iNummelm et all 



(200lb; 4 ref. lBergin e t al. (2005); 5 data collated bv lGibb et alT feoOO) 



One possible rea son for the failure of our models is a lack 
of gas-phase carbon: ISmith et al.1 (2004) tested increases in the 
C:0 ratio for the gas-phase model, and found improvements 
in the agreement with TMC-1CP, which by all accounts ap- 
pears to have an unusual chemistry. New estimates of carbon 
and oxygen abundances JSofia et al.lll994t iMeyer et al.lll998l) 
are both somewhat higher than those used here, but are un- 
lik ely to alter the C:Q r atio by a very great degree. As found 
bv lGwenlan et al. ( 2000l) . a much higher cosmic -ray ionisation 
rate may also contribute to the unusual chemistry of TMC-1. 
In any case, it is likely that we are lacking a crucial element in 
any attempt to model TMC-1 CP. 



4.3. Grain mantle composition 

The determination of ice composition in dark clouds is depen- 
dent on the absorption of infrared radiation from background 
field stars. Unfortunately, no data exist for lines of sight directly 
towards L134N or TMC-1 CP. However, a number of studies 
have concentrated on field stars behind the Taurus complex of 
dark clouds. Such data may give some indication of the likely 
grain mantle composition in the objects we have modelled. We 
compare our results with the ice composition determined to- 
wards background star Elias 16, which lies in the general vicin- 
ity of TMC-1 CP. Visual extinctio n along this line o f sight is 
determined to be around A y = 24 ( Whi ttet et al ] l2007l) . 

Table[5]shows the observed composition of icy mantles as a 
percentage of detected water ice, towards Elias 16. As the mod- 
elled grain mantle abundances are generally not very sensitive 
to the new desorption mechanism, we present results for model 
M2, chosen as representative of all models, at time 2.8 x 10 6 
yr, corresponding to the late-time peak in the agreement of gas 
phase species with observations of L134N, and at time 2 x 10 5 
yr, corresponding approximately with the early-time peak in 
agreement with both L134N and TMC-1CP observations (see 
figures |8] and [9). 

Whilst water abundance is fairly well-reproduced, there 
is clearly no strong agreement at either time-value for the 
strongest carbon-bearing constituents. CO and CO2 are under- 



represented in the model, with CH 4 comprising much of this 
carbon. Methanol is also somewhat more abundant in our mod- 
els than the upper limit of the observations suggest; however, 
a moderate reduction should have no major effect on the gas- 
phase abundance through the new mechanism. Formaldehyde 
is not detected, whereas quantities on the order of a percent are 
formed in our models. Ammonia is in acceptable agreement 
with observations, and H2O2 and OCS fall comfortably below 
their upper limits. 

The fact that the models do not well reproduce the main 
features of the observations is perhaps not surprising; previ- 
ous models using the canonical dark cloud physical condi- 
tions have not t y pically been very successful in this respect. 
Ruffle & Herbst ( 200 1 bl) attempted to rectify this, running a 
grid of models at various temperatures and densities and adopt- 
ing different grain surface binding energies and diffusion barri- 
ers. They found that greater temperatures and/or greater densi- 
ties could account for the particular surface composition seen 
in Elias 16, especially in reproducing CO2. We might there- 
fore suggest that the observed ice composition is representa- 
tive of a marginally higher temperature, say 15K. In our model 
this would lead to lower residence times for atomic hydrogen, 
resulting in less efficient hydrogenation. Our model adopts a 
set of binding energies and diffusion barriers representative of 
amorphous water ice; our atomic hydrogen diffusion barrier is 
a little higher than the lowest used by Ruffle & Herbst (225 K 
cf. 200 K), and our binding energy is rather higher than their 
lowest (450 K cf. 373 K). Our model also uses their low activa- 
tion energy (80 K) for the CO + O — > CO2 surface reaction. A 
lower hydrogenation rate (due to a higher temperature) would 
result in less CH 4 formation, and would allow CO2 formation 
to compete with hydrogenation, as in the Ruffle & Herbst mod- 
els. However, we might then argue that the Elias 16 observa- 
tions are not a fair comparison, as temperatures may well be 
different from those of the regions we are modelling. 

The apparent disagreement between the model and obser- 
vations should also be put into the context of our comparison 
with gas-phase abundances. Assuming a factor of 10 criterion 
for a successful match with observations, our ice results would 
appear much more acceptable, especially for earlier times of 
agreement. 

Nevertheless, the strong surface CH 4 abundance we find 
is not representative of observations of various other objects, 
typically reaching no more than a few percent of wat er ice 
dGibb et al.ll2000l) . As suggested by iBergin etail d2005l) . sur- 
face CH4 formation may be suppressed by a degree of gas- 
phase CO formation at earlier times, locking up carbon before 
significant grain mantles form. This would also favour greater 
CO2 formation, occurring in tandem with the hydrogenation of 
oxygen to H2O. Su ch a scenario agree s with the best fits to the 
ice observations of Ber gin et al.1 (|2005), which indicate that the 
majority of the CO2 resides in the deeper, polar layers. So the 
initial stages of evolution of the gas may be crucial to the ulti- 
mate composition of the ice, and indeed a static model such as 
this one may not be entirely adequate. 
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5. Conclusions 

The new desorptive mechanism, based on RRK theory 



(Hol brook et al.l 119961) . allows the chemical energy released 



from grain-surface addition reactions to break the molecule- 
surface bond of the product, with a small probability on the 
order of 1% governed by a parameter a. Under dark cloud 
conditions, each of the elements C, O, N and S (plus CO) 
and their various hydrides may participate in fast grain-surface 
hydrogen-addition reactions, which may therefore result in des- 
orption via the new mechanism. 

We delineate two main effects of the mechanism; an early- 
time and a late-time effect. At earlier times (t <~2 Myr), the 
appearances of gas-phase abundance profiles are in general 
not qualitatively altered, although in the case of a few select 
species, in particular, hydrogenated species such as methanol, 
the quantities are strongly enhanced. At late times, when de- 
pletion would otherwise be strong, the new mechanism acts to 
maintain gas-phase chemistry, by re-injecting the various hy- 
drogenated species - species that may achieve very large abun- 
dances on the grains. The re-injection occurs via photodisso- 
ciation of stable surface species into radicals followed by ad- 
dition reactions, which lead to desorption via the new mech- 
anism. This process feeds the formation of other species that 
were well-represented in the gas-phase up until the onset of 
strong depletion. Because of this effect, full-scale depletion is 
not attained, and a steady state is reached. It is important to 
note that in the case of CO, the strongest effect on both its 
gas-phase and grain-surface abundances, at late times, is not 
its own re-injection, but its gas-phase formation from surface 
hydrocarbons. Indeed, one of the most striking aspects of the 
late-time effect is the high abundance attained by grain-surface 
hydrocarbons. Our treatment of carbon chains in the gas-grain 
code, however, is certainly not comprehensive. Therefore, the 
abundance of these species may be overestimated, although we 
should expect some organic residue. However, the absence of 
observational evidence in dark clouds of solid hydrocarbons 
more complex than CH4 may suggest that the long timescales 
needed to acquire a significant hydrocarbon content in the grain 
mantles are not dynamically realistic. This would indicate that 
dark clouds, or the physical structures within them, are inher- 
ently short-lived, on the order of a few million years or less. 

Our new models show a good level of agreement with ob- 
served abundances for L134N. All of the models Ml - M3 
show improved agreement over MO, the model without the new 
mechanism. It is hard to constrain the value of the a-parameter 
by comparison with the observations of L134N; the net varia- 
tion in abundances between models Ml - M3 is not great. On 
the basis of methanol, and the general level of agreement with 
L134N, we would suggest a value of around a - 0.03 as the op- 
timum, and probably the maximum. (Such a value would sug- 
gest an optimum age for L134N of ~3 Myr, although we find 
that a much lesser age may also be plausible in this regime, 
perhaps as short as ~10 5 yr.) However, the most credible de- 
termination of a must ultimately be an ex perimental one. The 
simulations of Kroes & Anderssonl ( 20061) are encouraging in 
that they exhibit an effect of a similar strength for H2O. But 
due to the small size of the expected desorption fraction (~ 1% 



or less), the mere detection of the effect may be beyond exper- 
imental means at this time. 

Contrary to the case of L134N, observations of TMC-1CP 
are not well-matched by any of our models; the increasing 
agreement for continually higher a indicates that our models 
are not capable of producing a good match at any time, proba- 
bly because we lack some crucial information. A model that in- 
cludes dynamics/structure in the cloud may be required. Rather 
than a quiescent cloud, the TMC-1 ridge may be better repre- 
sented b y a superposition of clumps in varying stag es of evo- 
lution dGarrod et alj|2005[ l2006bl IPeng et alJI 1998b . A possi- 
ble explanation for the unique chemistry of TMC-1 is the ex- 
plosive desorption of grain mantles caused by MHD activity 
(iMarkwick etallbOOOh . 

Our models do not show strong agreement with the ob- 
served ice compositions determined towards Elias 16. In par- 
ticular, our models produce significant amounts of CH4, and 
much lower CO and CO? abundances. We suggest, after 
Ruffl e & Herb st (2001b), that a marginally greater temperature 
may better reproduce the measured ice composition; such may 
not be representative of the regions whose gas-phase abun- 
dances we model, however. Alternatively, greater gas-phase 
production of CO prior to signi ficant man tle formation may 
also improve agreement (Berg in et al 1l2005b . 



Qiang, Cuppen & Herbst (in prep.) have recently stud- 
ied cold cores with a new gas-grain approach, in which the 
grain chemistry is treated by a stochastic method, which is able 
to distinguish contributions from individual monolayers. With 
this approach, the observed abundances of CO2 and methanol 
ices in Elias 16 are fit reasonably at a time of 10 s yr whether 
one starts from fiat or rough olivine surfaces. As in earlier 
work, a temperature of 15 K results in improved agreement. 

Importantly, the inclusion of the new desorptive mecha- 
nism, over the range of a-values, does not drastically alter the 
observable chemistry of our models of quiescent dark cloud 
regions, but generally improves their agreement with observa- 
tions, while suggesting somewhat closer agreement at longer 
timescales. Whilst comparison with dark cloud abundances has 
been largely fruitless in strongly constraining the a-parameter, 
a model of pre-stellar cores may produce more success. These 
cold, dense regions are typefied by their strong CO depletion 
(e.g. iRedman et al. I l2002l) . Therefore, whilst for our quiescent 
dark cloud models the time at which depletion strongly takes 
hold is dynamically very late, pre-stellar cores must by defini- 
tion reach such a state within a dynamically achieveable period. 
We should compare chiefly against CO, due to its strong and 
well-documented depletion levels. A model that includes deu- 
terium chemistry would also be instructive due to the strong 
dependence of deuterium fractionation on the level of CO de- 
pletion. The degree to which the new mechanism may hold off 
depletion will be of interest in determining the ages of such 
objects. 
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